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Abstract 

In this paper we construct and approximate breathers in the DNLS 
model starting from the continuous hmit: such periodic solutions are 
obtained as perturbations of the ground state of the NLS model in 
H^iW^), with n — 1,2. In both the dimensions we recover the Sievers- 
Takeno (ST) and the Page (P) modes; furthermore, in also the 
two hybrid (H) modes are constructed. The proof is based on the 
interpolation of the lattice using the Finite Element Method (FEM). 

1 Introduction 

In this paper we study the problem of constructing breathers in the one and 
two dimensional discrete nonlinear Schrodinger (DNLS) equation starting 
from the continuous limit. 

The breathers we construct are critical points of the Hamiltonian func- 
tion constrained to the surface of constant i"^ norm. Such critical points 
are obtained by continuation from the continuous model constituted by the 
nonlinear Schrodinger (NLS) equation. The connection between the discrete 
and the continuous system is obtained by using the finite elements (FEM). 
This allows to identify the phase space of the discrete system with a subspace 
of the phase space of the continuous system. 

For example, consider the one dimensional case. The space of the finite 
elements is constructed as follows: first we associates to the j-th point of 
the discrete lattice a continuous piecewise linear function Sj{x), whose value 
is 1 at X = J and which vanishes for |x — j| > 1 (see Fig. [1]). To a sequence 
ipj, we associate the function Tp{x) := ipjSj{x/fi), where /u > is a small 
parameter representing the mesh of the lattice. The space generated by the 
functions Sj{x/ ^) will be denoted by 



Once this is done one can compare the functionals of the continuous 
system and those of the discrete one. In order to do this, denote by He and 
Nc the Hamiltonian and the norm of the continuous system, and consider 
the restriction of such functions to the space of the discrete system S^. By 
the standard theory of integration one can say that the restricted functionals 
are close to the Hamiltonian i/^ and the norm of the discrete system. 
So the idea is to consider a non-degenerate critical point of the functional 
of the continuous system, a critical point laying close to the manifold of the 
finite elements and to continue such a critical point to a critical points of 
the discrete functional. 

However there is a delicate point in the game: namely that the difference 
between the discrete functional and the continuous one should be small when 
the phase space is endowed with the energy norm. This turns out to be true 
thanks to a special property of the finite elements: the fact that one has 



with no error. Due to this property the difference between the continuous 
and the discrete functional turns out to be a functional which is small and 
smooth on the energy space. This allows to apply the implicit function 
theorem and to continue critical points of the continuous system to critical 
points of the discrete one. 

In order to be concrete we study in detail a one dimensional and a two 
dimensional model. We use known results on existence and non-degeneracy 
of the ground state of the continuous system in order to apply the above the- 
ory. In these paper we construct two (rcsp. four) kinds of discrete breathers 
in the l-(resp. 2) dimensional case, which are the continuation of the con- 
tinuous breather. In order to avoid problems related to the translational 
invariance of the continuous system we work here in spaces of reflection in- 
variant sequences. Thus the breathers we find for the discrete system are 
reflection invariant too. 

In dimension one, the breather of the first kind is centered at a lattice 
site and corresponds to the so called Sievers-Takeno mode (ST), while the 
breather of the second kind is centered in the middle of a cell of the lattice 
and corresponds to the so called Page mode (P). In dimension two, besides 
the ST and P modes, we have two other localized solutions, usually called 
hybrid (H) modes since they are centered in the middle of one of the two 
face of the cell. 

As far as we know, the result of the present paper is the first one in which 
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the continuous approximation is used in order to construct exact breathers 
of a lattice model. In dimension 1 the method of spatial dynamics also 
allows to construct and approximate breathers (see [Jam03]). However such 
a method is strictly one dimensional, while our method in principle applies 
to any dimension. Existence of breathers was also proved variationally in 
[Wei99] and in [AKKOl], but such methods do not allow to approximate the 
breathers and only allow to find one breather for each model. Breathers in 
DNLS have also been widely studied numerically (see for example [KRBOl, 
CJK+08,FW98]. 

The main advantage of our method is that it is quite flexible and allows 
to directly deduce informations on the shape of the breather starting from 
the continuous limit. 

We recall that the possibility of using the continuous limit in order to ap- 
proximate the dynamics of discrete systems has been widely investigated, in 
particular we recall the papers [BCP02,Sch98,KSM92,SW00,BP06,BCP09] 
in which an approximation valid for long but finite times and the papers 
[FP99, FP02, FP04a, FP04b, HW08, MP08] where an infinite time approxi- 
mation has been obtained. 

The plan of the paper is the following. In Section [2] we present the result 
and motivate our continuum limit approach. In Section [3] we formulate in 
Theorem 13. II the Implicit Function Theorem applied to our problem and in 
Section [4] we construct the FEM to interpolate the discrete model and we 
verify the hypothesis of Theorem 13.11 

2 Main result. 

We study here the discrete focusing nonlinear Schrodinger equation (DNLS) 
in M" with n = 1,2 

iV-i = -^(AiV')i - l^il^P^i , IGZ", \<p<- (1) 
/i^ In 

where Ai is the n-dimensional discrete Laplacian defined by 

(AiV')j := V-j+i + V'i-i - 2V'j, 
(AiV')j,fc := (V'i+i.fc + i^j~i,k - 2V'j,fc) + (V'j.fc+i + V'j.fc-i - 2V'j,fc), 

and /i is the lattice mesh. In particular we look for solutions of the form 

Ut) = e-'^Vi- (2) 
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Then the sequence ipi fulfils 



and thus it is a critical point of the Hamiltonian function 



1 IV'j-V'll 



1 



2p+2 



|j-i|=i ^ lez" 
constrained to a surface of constant value of the norm 

iV,:=^«^|V;,|2, 



(3) 



(4) 



(5) 



where the factors /u" have been inserted for future convenience. The main 
result of the present paper consists in showing that such a solution can 
be constructed and approximated starting from the continuous model 
constituted by the Nonlinear Schrodinger Equation (NLS), namely 



iip = —Alp ■ 



\2p 



(6) 



More precisely, consider the Hamiltonian and (the square of) the 1? 
norm N^^ given by 



|2p+2 



(7) 



then a periodic solution ^{x,t) = e tpix) of ^ fulfils the following con- 
tinuous approximation of ([3]) 

Xijj = -Alp - 



|2p 



(8) 



According to classical results on (HD (see [BL83, BLP81, CGM78]), there 
exists a unique real valued, positive, radially symmetric and exponentially 
decaying function ipc which realizes the minimum of Hc\nc=i- For example, 
in the case n = 1 and p = 1 it can be computed explicitly 
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ipdx) ■= -^sech ( — 



(9) 



If we interpret the discrete functional H^, as ^-perturbations of He, Nc 
and we restrict to a class of "even" functions in order to remove any possible 
degeneracy of the minimum ipc, then we can continue the solution ipc of ^ 
to a solution Tp{n) of ([3|). 

In order to state the precise result we are going to prove, we first need 
to define the configuration space for Tp\: 
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Definition 2.1. The space ^^(Z", 
with the norm 



will be denoted by when endowed 
1 



2_„(V,-AiV^)2,. 



(10) 



Theorem 2.1. For any fi small enough and ^ < p < ^ there exist 2" distinct 
real valued sequences V'i(^) ''^^^c/i are solutions of ([3|). Such solutions are 
even sequences ip^i = V'l lying on the surface = 1. One has 



\\lP' - 

where ^* is defined by 



Q <Cl^l, sup \^l;lit)-mt)\<C2^^ 



3 ri_ 

2 1 



fill 



IGZ" 









= i'c {m + f ) , 






\ \ 


:= V'c (w, AiA; + f ) , 








:= V'c (w + f + f ) 



n 



2. 



2.1 Comments. 

1. the first of ([TT]) is not empty since by using (jTTl) we get 

II^IIq.-!- 

Moreover, we stress that by its definition the approximating sequence 
is bounded uniformly in [i 



but is localizec|3 on an increasing interval [— /c, fe]" with A; ~ 1///. 
2. the first of (jlip immediatly implies 



\riit)-nit)\<c2fi' 



VI, 



an estimate which is empty in the case n = 2. Lemma 14.71 in Section 
4.3 is necessary to improve the above result. We do not know whether 
the exponent f — § is optimal or not. 



^We can fix the set where is localized as := {1 G Z" s.t. >I>j > i*o} 
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3. We stress that the problem ^ is equivalent to the /^-independent one 

Vi = -(Ai^)i - \ipi\^P^i (12) 

with the constrain 

This can be seen by the scaling 



and observing that 



2 2 

n > <^=^ p < —. 

p n 



3 The Implicit Function Theorem. 

The situation we will meet is summarized in the following abstract scheme. 
Let TChe a Hilbert space, and for any fi, let be a subspace of Ti.. Let He G 
C'^{TL) and Nc S C°°{TL) be two functionals, with Nc being a submersion. 
Correspondingly we define 

5 := {V G W : iVcW = 1} . 

Then we define the "discrete" objects: let //^ := H^^^^ G C^(<S^) and 
N(i := Nf:^,^ G C°°{£^) be functionals depending smoothly on two addi- 
tional parameters ei,e2- Define 

5,,,^ := {iPGSf, : iV.2,M(V') = 1} • 

We make some assumptions. 

i. There exists ipc GT~i which is a coercive minimum of Hc\q, namely it 
is a minimum and fulfills 

d^H,\^{i^,){h,h)>C\\h\\' , heT^^S; (13) 

moreover 

d{£^,ijc)<Cfi , (14) 

for all fi small enough. 
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Let ipo £ Sfj_ be such that \\ipc — ipo\\ < C/i and let U C £^ he an open 
neighborhood of ipo then we assume 



11. 



||-f^ei,/i - -f^o,At|lc2(w) - C'(ei + ^), Ho^fj. :— Hc\^^ (15) 
for some large enough k 



Theorem 3.1. Under the above assumptions, for any ei, e2, /i small enough, 
there exists a unique 1^^1,62, which is a coercive minimum of Hd\j^ ^j^. 
Moreover one has 

< + ei + ea) . (16) 
Proof. The result is local, so we restrict to a neighborhood of V'c- Define 

5o,^:={VG^:^ : iVo,^(V') = l} = cSnf^ (17) 
and take ipo G 5o,^. Remark that, due to smoothness of He one has 

\diH,\^Ji^Po)\\<Cfi. (18) 

By coercivity (|13p and Lax-Milgam Lemma, the second differential 

defines an isomorphism bounded together with its inverse uniformly with 
respect to all the parameters. 

From assumption ii, there exists a local isomorphism 

which satisfies 

\\Ie2,t.-Id\\ck < C{e2 + ^i) . 

The statement is then equivalent to the existence of a coercive minimum of 
-fTg^^^ o 1^2,^1- To get it remark that 

oT,,,^ = 774^^^ + 0(ei + £2 + /i) . (19) 

Due to (fTSll and (fTSll . the Implicit Function Theorem applies and gives the 
result. □ 
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4 Applications to breathers 

In order to avoid gauge and the translational invariance of the problem, in 
particular of the continuous system, we will work in a space of real valued 
functions "invariant" under the involution 

S : ip{z) ^ i^{-z), z G M". (20) 

More precisely, in ff-'^(M^,M) we will consider functions fulfilling 

f |^(^)_^(_^)|2 = ||^_5^||2^=o, (21) 

which is equivalent to (|20p almost everywhere and is a condition well defined 
in Fi(M^R). 

Lemma 4.1. Let tpc be a solution of ([SD with p < ^ and let 

W = {V G H\W'\R) : Wi; - Sij\\l, = o} =: H^, 
then assumption (jl3p of Theorem \3.1\ holds. 

Proof. This Lemma directly follows from Proposition D.l of [FGJS04] by 
remarking that T^^S C X, with X defined in the statement of Prop. D.l. 

□. 



Remark 4.1. We stress that the constrain (|2ip is "natural" for the problem 
([3|), since (j20p is a symmetry for both the Hamiltonian and the Norm 
([5]). Hence, a critical point for the restricted problem is also a critical point 
for the original problem. 

In the following subsections we construct the linear manifold of the 
finite elements, and prove the estimates (jl4p and (jlSp for the two considered 
applications. We deal with the ST-breather, since the other ones follow by 
small changes in the definition of £f^. 

4.1 The case n = 1 

Let 1 = j and define the sequence of functions Sj{x) by 



Sj{x) 



0, if \x - j| > 1 

x-j + 1, if j-l<x<j (22) 
^ —X + j + 1, if j < X < j + 1 
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j-1 j j+1 



Figure 1: Finite element basis Sj{x) and interpolating function ijj{x). 

and, to a sequence ipj G Q^, we associate a function 

^(x) :=^Vi5i(^//")- (23) 
3 

On the interval Tj := + 1)) the above function reads 

^(x) = (x-wO^^^^^^^^ + V-i- (24) 

Definition 4.1. M^e denote by £^ the linear space composed by the functions 
of the form ([23]) with ipj G Q^. 

The following Lemma gives the equivalence between the function space 
£^ and the sequence space Q^. 

Lemma 4.2. Let ^' G then 



Moreover 



[ ^l = ^^{i;^^,-i^^Y. (25) 

l^Y.^]= f^' + ^f ^Idx. (26) 



Proof. Let us first decompose M = Uj^zTj- The weak derivative of ^' is 

^.(x) = ^ ^^^-^^-^^-) xto,M.+i))(^) 
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which gives immediately 



nix) 



■XtAx) 



E 



^^^±i-r^XT,(.), (27) 



since XTj{x)xTi{x) = for i / j. From J^XTj{x)dx = fi one gets 



(28) 



If we plug (p^ in the integral ||^||^2, a direct computation gives the estimate 
(I26D. □ 



Proposition 4.1. Lei ^' G f ^ he as in (j23|l an(i /ei us define 



g+2 
'J I ' 



if q > ^ then Rq G C^iSfj,) and for any bounded open setU C <5^, i/iere exists 
CiJA) such that 

WRaWc^ipi) < Cfi. (29) 

Proof. The term Rq can be represented through the Euler-MacLaurin for- 
mula 

E/(J')= / f(y')dy+ [ fy{y)Pi{y)dy, Pi{s) = s - [s] - 1- (30) 
Indeed, if we set f{y) = |^'(/iy)|'?"^^, we have 



j;|vI/(;.i)|'?+2^^^/(j)=;. f \^{^,y)\'^+^dy + 



fi^ / ^{f,y)\^{fiyW^,if,y)Pi{y)dy = 

[ \^ix)\'^+^dx + fi [ ^{x)\^ix)\'^^^{x)Pi{x/n)dx. 
Jr Jr 



10 



Figure 2: Triangulation and finite element in the bidimensional lattice. 
Hence 

Jr 

A direct computation of the firsy and second differential shows that 

dRffi'^'m = n [ {^m''),hPiix/i2)dx + f, [ ^\^\''h,Pi{x/fi)dx, 

d^RH{'^)[h,h] = [ m'')xh^Pi{x/fi)dx + [ m'^{h^),Piix/f,)dx. 
Jr Jr 

The smallness is represented by the prefactor fi: so Sobolev embedding 
Theorems and |Pi(x//i)| < 1/2 yield 1^. □ 

We have thus verified the assumptions of Theorem 13.11 which implies the 
existence and the estimate of the ST- mode for the case n = 1. The same 
statement for the P-mode follows by a translation of the basis of i?^ 

^{x) := ViSj(a://U + 1/2), V'-i+i = V'j, (31) 

with Sj defined in (j22p . 

4.2 The case n = 2 

Let us take ipj^k G Q^- For each multindex 1 = (j, k), let us consider the func- 
tion Sj^k{x, y) which represents the exagonal pyramid of height one centered 
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in (j, k) whose support is the union of the six triangles of figure [2j More pre- 
cisely we define Tj~j^ the triangle whose vertexes are {j, k), (j + 1, A;), (j, A; + 1) 
and T~f, the one whose vertexes are {j, k),{j — l,k), (j, k — 1). Hence, for 
example, on T^^^ the function sj^k represents the plane in 

Sj,k{x, y) = -X - y + j + k + 1. 

The set of functions {sj^k{x/fJ',y/fJ')}(j^k)ei.'^ is a basis which generates a 
piecewise linear function interpolating ipj^k 

*(a^>y) := X] '^j,kSj,k{x/fi,y/fi). (32) 

Notice that on the triangle T^j^ the function ^' is the plane 

^{x, y) = ipj^k + {x- fij)^^ + {y- fik)^y. (33) 

Definition 4.2. We denote by£^ the linear space composed by the functions 
of the form ([32]) with ipj^k ^ Qfi- 

The following Lemma gives the equivalence between the function space 
£^ C and the sequence space Q^. 

Lemma 4.3. Let ^ €z £^ then it holds true 



(34) 



(j,fc)ez2 -"^ {j,k)&:^ 

Moreover 



^ 7lR2 6 ' 



(35) 



Proof, from (j33p we have that on each triangle T^^ it holds 

Formula ([35]) follows from a direct computation as in Lemma |4.2[ □ 



The next three Lemmas provide the proof of the following main 
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Proposition 4.2. Let € £^ be as in ([32]) and let us define 
G,{^) := [ \^\'^+^ 

if Q ^ ^ then Rq G C'^{£^) and in any open set U C £^ one has 



Lemma 4.4. Under the assumptions of Proposition \4-^ one has 

\RGm<Cf,\m'H,\ 

Proof. Let us set f{x,y) = \^ (x , y)\'^^'^ and let us take (x,?/) G Tj^^, then 
we can use a Taylor expansion with integral remainder 

fix,y) = f{fJ-j,fik) + {x-fij) / fx{'y(t,x,y))dt+{y-nk) / fy{'y(t,x,y))dt, 

Jo Jo 

(36) 

where 

7(t, X, y) = (tx + (1 - + (1 - t)/iA;) t G [0, 1] 

is the segment connecting {p,j,fj,k) with and lies in the triangle T^f^. 
Hence 

= ^f{nj,nk) + 

+ / {x - nj) fx{'yit,x,y))dtdxdy + (37) 

+ Ay-fJ-k) fy{j{t,x,y))dtdxdy. (38) 



By the initial definition of f{x, y) one has 







-1 IV'iii.fc - ipj. 



\d.f\ = (g + 2)|^r+^|M/.| = (g + 2)|^&r^ '^^^^'^ , (39) 

= (g + 2)|vI/|«+i|M.,| = (g + 2)|*r+iM^^^±i^^. 

(40) 
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Since 

,2 



fix,y)-^fij,k) 



we can use (j37|) and ([38]) to estimate 



fJ- ^ \fx{l{t,x,y))\ + \fy{j{t,x,y))\dtdxdy. 



< 



By inserting (j39p and observing that 

o 7 = ^'(7(t, X, y)) = Vi,fc + t(x - + i(y - (41) 

one may remove the integration along the segment 

where we set Sj^ki^, u) '■= {x — fJ-j, y — fJ-k) and 

(V^, 6j^k) = {x- //j)^'^. + {y- fik)^y. 
Notice that from ()33p one has 

\i^j,k\ < \^{x,y)\+ij{\^^\ + \^y\) 
thus it is possible to estimate the argument of the integral as follows H 
+ |(V*,<5,-fc)|]^+i < 4^ b(x,y)|^+i + ^'^+1 (1^,1 + \^y\y+' 



^We here remind the inequahties for a,b > 

{a + by <2'{a' + b''), 0<s<l, 

and 

{a + by <2'-\a' + b''), s>l; 

the second follows easily from the convexity of the function g{x) — x", x £ R^. The first 
is a direct consequence: 

« + by = \ < 2" + < 2" a" + b"). 

{a + b) - \{a + b) {a + b)J - ^ ' 
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j,k& 



hence recalling also that l^'^^l + l^'yl < \/2|V^| one has 

= MC(g) X / |VvI/|[|v&(x,y)|''+i + ^9+i|V^r+'] < 

/ |V*||^'(x,y)|5+^ + ^'?+l / 

The first sum gives 

< Ci||V^I/||^.||^||^+\ 
Using Lemma 14.31 the second instead gives 



19+2 

I V :i: 

n21 1+f 



j,fcGZ2 

< c,^,\\vnt'■ 

Collecting the above estimates we obtain 
which finally gives 



/ |^I/(x,y)r+2_^2 ^ |^^.,|'^+2|<^c 11*11 



g+2 
HI 



j,fcez2 

□ 
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Lemma 4.5. Under the assumptions of Proposition one has also 

II^GW|li(,,,M)<^/^ll^llH^ (42) 
Proof. A direct computation easily gives for any h ^ £^ 

with obviously hj^k = h{iJ,j,fik). In order to estimate 

\\G,{^} - - sup 

we need to control 

\G',mh]-G',mh]\ = \ [ (43) 

We proceed in the same way as in Lemma 14.4^ exploiting the fact that also 
h{x,y) £ We thus define in this case f{x,y) = y), 
so that 

i/xi + \fy\ < q\^\'{\^.\ + i^yim + i^r+i|v/i|. 



We recall that 



\hoj\ < +2/i|V/i|, 



hence (I43p can be split into four terms 

+ [ i^-i^+^iv/ii + 
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1. using Schwarz and (V|^'|'?+^)2 = V^'Pl^'p'? we get 



'\\h\\H^< 



2. using Schwarz and ^ s > 2 



< C2\\vnit'\\h\\j,.; 



3. as in 1. 



1 ii^ii^i^ ; 



4. as in 2. 



< Ci llV^ll^f ||/.||^. . 

Collecting we get 

\R'Gmh]\<c{qh\m'j^.'\\h\\^,, 

hence the thesis. □ 

Lemma 4.6. Under the assumptions of Propositions one has also 

Km\L\e,,^)<Cf^\m'H^- (44) 
Proof. Also in this case, a direct computation easily gives for any h £ £^ 

GUm,h]= [ \^\'h^ G'm[h,h]=f,'Y.\i;,,k\'hl,. 
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In order to estimate 

we need to control 

\G'^{m,h]-G'^mh,h]\ = \ [ \^\''h'-n'Y.\^,,,\''hl,\. (45) 

We proceed as in the previous Lemmas, by defining 

f{x,y) := \^{x,y)\ih?{x,y), 

so that 

\fx\ + i/j/i < + 2mhvh\. 

We distinguish the case g = 3, which is easier, and g > 3. 
g = 1 In this case we have 

E/+ ri/-°7l + l/,o7| < 2/ \Vnh? + 

+ V / |v^||v/i|^ + 

+ / |^'||V/t|(|/i|+/x|V/i|) + 
+ /X / |V*||V/i|(|/i| +Ai|V/i|). 
The thesis can be obtained using Schwarz and observing that 

JM2 M . , A* 



4 

L2 



g > 1 The steps are the same as usual; the only difference is that we have 
to deal with 

I^P'^IV^'P, <T>0, 



but it is enough to notice again that the integral above is the (square) 
norm of Vl^'|^+'^, thus 



7m2 
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□ 



This concludes the case related to the construction and approximation 
of the ST-mode. The other three modes (the P-mode and the two H-mode) 
are obtained by translation of the basis Sj either in one or in both the two 
directions. 



4.3 Proof of Theorem \27\\ 

We begin with the following 

Definition 4.3. Let n = 1,2 and consider ip G ^ on 8^. We 

define 

I^^■.^l:^I^^^=Y,^{^J^^)s\{x/^J), xew^ (46) 



the projection of H (W^) on By classical results on polynomial approxi- 
mation in Sobolev spaces ( Chapter 4 of [BS08]) one has 

m^.^-nm<c^iu\\H^- (47) 

We need also a simple lemma to obtain the second estimate of (|lip 
Lemma 4.7. For any 1 G Z" we have 

1 n 

1^,1 <2/.2-2 llV^llg (48) 



Proof. We write the proof for the case n = 2. The case n = 1 is simpler. 
Denote 1 = (j, k); one has 

h h 

m= 

which gives 



m=— oo 



m=— oo 



sup < 4 / ^ ^jI^ j^ / (V'm+l,fc - ^m,kY < 



mG2 



mm 



< 



i2H 



2 4 



1^' 
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and 



□ 



Now we easily verify the hipothesis of the abstract Theorem 13. 1[ First, 
we define ipc as the (smooth) solution of dH]) and ipQ = i^H^ipc, with v such 
that A'^o^^(^o) = 1- Then, condition (I13p follows from Lemma 14.11 while 
condition (jl4p comes from the above (j47p . Finally, requirement ii is given 
by Lemmas 14.21 and 14.31 and by Propositions 14. l l and 14.21 This directly gives 
the first of pT|) . The second of ([IT]) is a byproduct of either the first and 
Lemma [121 indeed 

□ 
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